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ABSTRACT 


This article deals with the application of a new 
graphical system, known as Graphical Rational Patterns (GRP), to 
geostatistical data. The first portion of the paper provides a 
brief review of the system for readers who are unfamiliar with 
it. This is followed by a series of annotated illustrations based 
on recent Canadian data. The paper emphasizes that the GRP system 
is particularly suitable for producing large numbers of inexpen- 
sive drafts or experimental maps by non-cartographers; as such, 
the system contributes to geostatistical analysis undertaken by 
demographers, statisticians and others. An attempt is made to 
substantiate that the GRP system (7) introduces improved arrays of 
sympols (e.g,,with regard to homogeneity and "rationality"); (7%) 
provides graphs and maps with improved "readability" and transmit- 
tance of general trends; (iit) suggests solutions where conven- 


tional methods are inapplicabde or ineffective. 


1. INTRODUCTION 


One hardly needs to convince statisticians that in 
analysing statistical data "one picture - that is, a chart, 
graph or map - is worth a thousand words". This applies particu- 


larly to geostatistics? where the difficulty of conceiving 
I. The term "geostatistice", not yet used universally, refers to the analysis 
of geographical-statistical dietributione. 
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trends embodied in statistical tables is compounded by the need 
to envisage geographical contexts as well. In this connection, 
conventional? cartographic techniques are inadequate for treat- 
ing several important categories of series, e.g., double entry 
tables showing cross-classifications>. Even more disturbing, in 
the author's opinion, i8 thé Tack of 4 géherdl tdol enabling 
analysts who are not professional cartographers to produce large 
numbers of inexpensive and easily prepared graphs for expe- 


riments, drafts, interim reports etc. 


The system known as Graphical Rational Patterns (GRP) 
provides some answers to these problems as well as to other defi- 
ciencies of conventional techniques, to be discussed later. How- 
ever, though the GRP was introduced several years ago by its in- 
ventor, R.Bachi(1968), the system has been applied only rarely, 
the official publications of Israeli’ statistics being an excep- 
tion [see for example the summary volumes 42 and 43 of the 1961 
census of Israel]; textbooks on cartography, with the exception 
of the one by Bertin (1967,p.392),seem to ignore the GRP comple- 
tely. 


The aim of this paper is to draw the attention of Cana- 
dian statisticians and demographers to the potential of the GRP 
in data analysis. To this end, a series of annotated charts and 
maps illustrating selected features of GRP presentations, are 
shown in Section 2. For the benefit of readers who are not 
acquainted with the GRP system the remainder of this section 
reviews its basic properties. A detailed discussion may be found 
in Bachi(1968) and Yam(1972). 


Essentially, the GRP system consists of a series of 
graphical schedules suggested in lieu of conventional ones. The 
basic GRP schedule, termed the single GRP scale, provides 101 
graphical symbols which correspond to the integers 0-100 
(Fig.1). A GRP unit, that is, a black square of given area, 
represents the integer 1. Two, three, ... , nine GRP 


2. By "conventional" or "traditional" we mean the graphical techniquee one 
te Likely to encounter in text-booke, offictal publications, etc. 


3. Stereagrans (also known as "atepped etatistical eurfacea") may be considered 
an exception but they are quite impractical on the whole. See, for example, 
Zelinski (1966, p.8). 
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units are clustered in distinct configurations to represent the 
corresponding integers 2 to 9. The particular cluster configura- 
tions are selected to assist in identifying the given symbol, 
relative to others in the array. The integer 10 is represented 
by one black square, the area of which equals 10 GRP units, that 
is, ten times the area of the symbol corresponding to the integer 
1. The nine configurations employed to show the numbers 1 to 9 
are now repeated with the larger squares, thus forming the sym- 
bols for 10, 20,...,90 (compare the symbols for 5 and 50, as an 
example: the configuration is identical but the size of each 
"pbuilding block" is different). A combination of the units and 
tens thus yields a scale from 1 to 99. Zero is represented by a 
blank GRAF and the integer 100 by a solid black square, the area 
of which equals 100 GRP units. The symbols described above are 
subsequently referred to as GRP's or ‘GRP frames". The correspond- 
ing integers are termed "basic values". 


The basic property of the GRP scale is its rationality, 
in the following sense. Let b,, denote the proportion of the inked 
surface of a given GRP symbol (this is termed the "blackness" 
of the symbol); let n be the corresponding basic value. From the 
construction of the GRP scale it follows that: 


(1) by: be: ss) are be: Paes bE 1 OS) placate, «ag Lome oie 


A series of extensions and variations have been devised on the 

basis of the said single GRP scale. The main additional sched- 
ules include: (7) an extension of the 0 - 100 single GRP's to 

1000; (it) repeated linear GRP's, scale representing the integers 
0 to 20 (applied to presentation of traffic lines, for example) ; 
(iit) repeated areal GRP's, scale representing the integers 0 to 
tO (Fig.3); (iv) multiple GRP's for a simultaneous representation 


of two variables, each in the range 1 to 10. 


Additionally, the single GRP’s may be enclosed in 
frames of various shapes - square, round, oval, etc.- to 
introduce an additional qualitative variable or to indicate neg- 
ative numbers (Fig.6). Various colour schedules increase the versa- 


tility of the system further. F 
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1 — COMPLETE SCHEDULE OF SINGLE GRP®*Sy 


FIG. 


VALUES FROM 1 TO 100 


91 


81 


vel. 


61 


51 


41 


30 
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ILLUSTRATIONS“ 


2. 


2.1. MAPS OF POPULATION DISTRIBUTION - A SIMPLE APPLICATION OF 
SINGLE GRP'S 


According to a common traditional technique, maps show- 


ing the regional distribution of a papulation are prepared with 


the area of which is proportional to the 


circular symbols, 


The GRP system proposes to sub- 
4, This footnote ts given on page 236. 


population in a given region. 
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stitute single GAP’s for the circular symbols. A map of this 
kind is shown in Fig.2, where each GRP unit represents 100,000 
people in any given province of Canada, 1971. The reader in- 
terested in a general picture needs only notice the larger black 
squares representing 10 GRP units, that is, 1,000,000 inhabitants 
Thus the GRP map conveys an effective visual impression of the 
Major trends latent in the series. One advantage of Fig.2 over 
the traditional equivalent stems from the easy readability of 
GRP graphs: the simple correspondence between the GRP symbols 
and the populations represented assures the reader that decoding 
the symbols may be achieved easily, should he want to exploit 
this option. Experience has shown that with minimal practice, 
decoding may be exercised with little reference to the "key" to 
the map. Furthermore, since the GRP is essentially a language, 
experience acquired from one graph is applicable to GRP graphs 
generally. Fig 2 also illustrates how the GRP system facilitates 
the preparation of geostatistical maps. Firstly, planning the 
map is an easy task thanks to the large number of patterns in 
the array; secondly, applying the GRP's in the preparation 
proper is a matter of minutes. 


2.2. DENSITY MAPS - A SIMPLE APPLICATION OF REPEATED AREAL GRP'S 


GRP density maps substitute a homogeneous” setiof 51 
repeated areal GRP symbols for the conventional array of pat- 
terns - usually a heterogeneous one, consisting of 5 to 7 pat- 
terns. Fig.3 demonstrates a GRP density map by County/Census 
Division for part of eastern Canada, based on the 1971 Census 
(the raw data, for comparison, are reported in Table 1). In the 
author's opinion, the map confirms the claims made previously, 
‘concerning transmittance of trends, easy readability when 
desired, and assistance in map preparation by non-professional 


4. The illustrations shown in Section 2 were prepared with virtually no 
assistance of profeastonal cartographers. Only standard base maps, com- 
mxcially avatlable GRP's, conventional photocopying machines and computer 
print-oute were ueed, Thie method was adopted expresely to demonstrate 
the type of tentative graph or map that analysts may prepare rapidly in 
the course of their work. The graphs as shown should not be considered, 
therefore, as products suitable for final publication in research 
reports. 


5. "Homogeneous", in the sense that the GRP symbole are all constructed 
from the sane building blocks, t.e., the GRP unite. 
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2.3. GRAPHS REPRESENTING DOUBLE ENTRY TABLES 


While the previous two examples dealt with GRP maps, the 


objective of the present illustration is to demonstrate a general 


application of the single GRP's to statistical charting. Thevir- | 


tually uniform size of the single GHP frames renders them partic- 
ularly suitable for graphical representations of double-entry 


tables. As an example, consider the Canadian 


internal migration 


matrix by Region, as per Table 2 and the corresponding Fig.4. 
Though the contribution of graphs to analysis of such data is 


best appreciated when large matrices are dealt with e.g., internal 


migration in Ontario by County), the 5x5 matrix employed here is 


sufficient as a general indication. 


TABLE 1 — POPULATION DENSITY (PERSONS PER SQUARE MILE) BY 


CENSUS DIVISION (C.0.) OR COUNTY, 
P.E-1. AND NOVA SCOTIA, 1971 


Cre) eee UL Ry tO RIN Sli tny, C.D. OR DENSITY 
COUNTY COUNTY 
NFLD. 366 NeSe 38.7 
1 60.6 1 17.7 
2 liel 2 29.8 
5 3.0 3 W352 
4 928 4 27.0 
3 12.6 5 21.2 
6 5.7 6 21.3 
7 10.5 u 726 
8 13.6 8 121.8 
9 4.1 2) 24.5 
19 0.3 10 14.3 
Lt 53.4 
12 34.6 
13 43.1 
PET. Divers 14 14.2 
------ 15 26.8 
1 28.5 16 18.3 
2 bDrek id 7.3 
3 66.0 18 30.9 


NEWFOUNDLAND, 


NEWFOUNDLAND 


SOURCE - STATISTICS CANADA, CENSUS DIVISION, 
(UNPUBLISHED DATA). 
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FIG. 


SOURCE? 


NEWFOUNDLAND (EXCLUDING LABRADOR)+ 


TABLE 1. 
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AND NOVA SCOTIA, 


3 -— POPULATION DENSITY (PERSONS PER SQUARE MILE) BY COUNTY OR CENSUS DIVISION, 
P.E.I. 


1971 


KEY: 
1 GRP UNIT = 10 PERSONS/SQUARE MILE. 
THE COMPLETE SCHEDULE IS AS FOLLOWS: 
==--SYMBOL-==— | —---—= VARI ABLE----—- 
PATTERN BASIC CLASS RANGE CENTRAL 
VALUE VALUE 

0 On0!— 409 (0) 
Be eee 1 5.0 - 14.9 10 
Rana 2 15.0 — 24.9 20 

3 25.0 = 34.9 30 

4 3520 - 44.9 40 
mse SELSE 5 45.0 - 54.9 50 
238 32% 

6 55.0 - 64.9 60 
a 7 6520 — 14.9 70 

8 75.0 - 84.9 80 
aoe 9 85.0 - 94.9 90 
[an 10 95.0 + (100) 
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Examining an internal migration matrix, analysts would typically 
transform the raw data (Table 2A) into percentages (Table 2B) 

and proceed to compute "expected" values on the assumption of 
“independence of marginal totals" (Table 2C). If m,, and é,; denote 
the "observed" and "expected" streansfrom i to j (in percentage 
form), respectively, then: 


2 we = Chom. #8 fs 
(2) 2 at ; he /100 


The next step may involve the computation of the ratios: 
“observed"/"expected" (that is, MAE, 5) 0 to indicate the degree 
of preference of various streams (Table 2D). Usually, the 
extreme cases (that is, the most preferred streams as well as 


the least preferred ones) attract the analyst's attention most. 


Figs.4A to 4C illustrate how these data may be presented 
in graphical form by using the single GRP scales. Referring to 
the main diagonal of Fig.4A, for example, it is readily observed 
that migration within Regions (according to the definition of 
Table 2) accounts for the lion's share of internal migration in 
Canada; from the diagonal of Fig.¢C it is learnt further that 
the "observed" streams are at least 2.5 times as large as the 
"expected" ones for all intra-regional migrations. This example 
suggests that GRP charts may assist the analyst in visualizing 
the trends latent in complex statistical tables. As mentioned 
in the introduction, no practical alternative is proposed by 


conventional cartography. 


2.4. MAPS OF SERIES SHOWING "REGIONAL AVERAGE INTENSITY 
RELATIVE TO A BASE POPULATION" ("INTENSITY SERIES") 


Let P; denote the population of region 7 as enumerated 
on one census date; let P; correspond to the population at a 
later census date, so that the percentage change A, is defined 


(3) A. = 100 (P} -P,)/P ; 


For graphical presentations of typical "intensity series", such 
as the one formed by the values Aye conventional choropleth 

maps are commonly proposed. This type of map and one of several 
GRP alternatives are demonstrated in Figs.5 and 6, respectively. 
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FIG. 4 ~ MATRIX OF INTERNAL MIGRATION, CANADA, 1961 CENSUS DATA 


A. OBSERVED STREAMS ( PERCENTAGES ) B. EXPECTED STREAMS ( PERCENTAGES } 


pe 
TO REGION - | 
a 
FROM ATLNT-| P.@. | ONT. | PRAIR.| B.C. | CANADA BECK ATLNT.| P.Q. 


REGION i a) EGION | 


e dl a ATLNT. . " 4 “ : =] 
cnee8 pads fe Fl 


— + } = + 
|| a = ie if 
] a a 
ONT. . “ i] . . 5 ry Z “oe a a is 


TO REGION 


| 
7 


ATLUNT. 


+————| 
PRAIR. ~ ae * at PRAIR« ~ o a 5 s mo 
at L mam + — r a 4) = FI 
B.C. . s a “a BeCe . 4 a ” . “a 
LS eS Se —— a —— we =" 
g B a * be 5 a a 
CANADA a ny CANADA * = 
= | o's a = = = | ee a — 
==: [ee ut ——J —— md 


C. RATIOS: OBSERVED / EXPECTED ( PERCENTAGES ) 


FROM 
REGION 
ATLNT. 


TO REGION 


CANADA 
=) 


KEY TO FIGS. 4Ay 4B AND 4C: 
1 GRP UNIT = 1%, EMPLOYING THE SINGLE GRP SCALE 
EXTENDED TO 300. 


EXAMPLES: 


CANADA 


249% 300E+ 


SOURCE: TABLE 2. 
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The major difference between the two maps is the vi- 
sualimpression generated by Labrador, and the question arises as 
to which of the maps transmits correctly the message latent in 
the data. Clearly, both maps attempt to represent the variable 
4 by means of the basic property of the graphical symbols, that 
is, by their "blackness" (as defined in Section 1). If b, denotes 
the blackness of the symbol representing Aes it is reasonable 
to require that the visual impression generated by region 7 
depend on b;- In the conventional choropleth map, however, the 
visual impression is created by weighting the b's by the map-are 
a, Of region zi. It is suggested here that (zi) the weighting fac- 
tors a; are irrelevant to the presentation of A; and hence the 
conventional practice is unjustified; (ii) the overwhelming vis- 
ual impression associated with regions such as Labrador (Fig.5) 
is only in part due to the variable A and is mainly a consequence 
of the region's area - hence the conventional practice may mis- 
lead the reader. The simplest GRP alternative applies one frame 
of the single GRP schedule per region, such that its b, is propor 
tional to A. No weighting factor is introduced and hence the 
area of the region does not affect the visual impression. An 
alternative GRP technique suggests the preparation of weighted 
Maps for intensity series. Accordingly, D; frames, rathen than a 
Single one, are assigned to region 7: unlike the choropleth maps, 
however, the weight P; is determined so that Baw Ps. As a 
result, and recalling expression (3), the total blackness in 
region t is given by 


« « fii 
(4) pb, = Pb, = PL-P.. 


Expression (4) implies that the overall visual impression gen- 
erated in region t is proportional to the magnitude of the inter 
censal population change (in absolute numbers). 


Reviewing Fig.6, which represents a weighted GRP map, 
one may easily determine that the value of 4 is higher in Labrdor 
than in any other Census Division, yet thé visual impression 
created by this region is not as overpowering as in Fig.5. On the 
other hand, and in contrast to Fig.5, the map leaves little doubt 
that the absolute change in Labrador is smaller than the one 
which took place in the Census Division of St.John's. Sum- 
marizing the foregoing discussion it is suggested that GRP maps 
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for intensity series are preferable to conventional presentations. 


3. CONCLUSION 


In the course of commenting on the four GRP maps 
presented above, an attempt was made to substantiate the proposi- 
tions that the GRP system (i) answers the basic requirement of 
any graphical system, namely, that the essential properties of 
the data represented be conveyed visually in an effective manner; 
(tt) offers an additional option of easy readability; (iit) 
provides new solutions in certain cases where conventional methais 
are impractical or otherwise deficient; (iv) enjoys additional 
advantages, theoretical (e.g., rationality) as well as aesthetic 
(e.g., homogeneity). Significantly, the practicality of the GRP 
as a system which may be applied by non-professionals was em- 
phasized. It seems appropriate, then,to conclude on a practical 
note, suggesting that statisticians and demographers experiment 
with this new system in the course of their work. 


RESUME 


Nous considérons l'application d'un nouveau systeéme graphique, 
connu sous le nom de Graphical Rational Patterns (GRP), a des 
données géostatistiques. L'article fournit d'abord une bréve 
révision du systeme et de sés principes. Ceci est suivi par une 
série d'illustrations cartographiques et par des commentaires. 
Toutes les cartes sauf une sont basées sur des données 
canadiennes (principalement du recensement de 1971), l'exception 
se rapportant aux Etats-Unis. Pour permettre au lecteur de se 
rendre compte des avantages du systéme GRP, des tableaux des 
données originales accompagnent chaque carte, et dans trois cas 
les cartes conventionnelles équivalentes sont ajoutées. L'artick 
accentue le fait que le systéme GRP est particuliérement adapté 
a la production d'un grand nombre de dessins peu couteux ou de 
cartes expérimentales par des personnes qui ne sont pas des 
cartographes; ainsi ce systéme contribue a l'analyse géostatis- 
tique entreprise par des démographes, des statisticiens, etc. 
D'autres caractéristiques du GRP comprennent: (i) l'introductim 
d'un ensemble amélioré de symboles ayant trait a l'homogénéite, 
“la rationalité", et le nombre de configurations; (ii) la 


préparation de cartes avec une lisibilité accrue et la mise en 
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evidence des tendances; (ttt) des possibilités de solutions ot les 
methodes conventionnelles sont inapplicables ou inefficaces. 
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